Conformational Space of the Translocation Domain of Botulinum Toxin: Atomistic Modeling and Mesoscopic Description of the Coiled-Coil Helix Bundle.

The toxicity of botulinum multi-domain neurotoxins (BoNTs) arises from a sequence of molecular events, in which the translocation of the catalytic domain through the membrane of a neurotransmitter vesicle plays a key role. A recent structural study of the translocation domain of BoNTs suggests that the interaction with the membrane is driven by the transition of an α helical switch towards a ß hairpin. Atomistic simulations in conjunction with the mesoscopic Twister model are used to investigate the consequences of this proposition for the toxin-membrane interaction. The conformational mobilities of the domain, as well as the effect of the membrane, implicitly examined by comparing water and water-ethanol solvents, lead to the conclusion that the transition of the switch modifies the internal dynamics and the effect of membrane hydrophobicity on the whole protein. The central two α helices, helix 1 and helix 2, forming two coiled-coil motifs, are analyzed using the Twister model, in which the initial deformation of the membrane by the protein is caused by the presence of local torques arising from asymmetric positions of hydrophobic residues. Different torque distributions are observed depending on the switch conformations and permit an origin for the mechanism opening the membrane to be proposed.

In Silico Conformational Features of Botulinum Toxins A1 and E1 According to Intraluminal Acidification.

Although botulinum neurotoxins (BoNTs) are among the most toxic compounds found in nature, their molecular mechanism of action is far from being elucidated. A key event is the conformational transition due to acidification of the interior of synaptic vesicles, leading to translocation of the BoNT catalytic domain into the neuronal cytosol. To investigate these conformational variations, homology modeling and atomistic simulations are combined to explore the internal dynamics of the sub-types BoNT/A1 (the most-used sub-type in medical applications) and BoNT/E1 (the most kinetically efficient sub-type). This first simulation study of di-chain BoNTs in closed and open states considers the effects of both neutral and acidic pH. The conformational mobility is driven by domain displacements of the ganglioside-binding site in the receptor binding domain, the translocation domain (HCNT) switch, and the belt α-helix, which present multiple conformations, depending on the primary sequence and the pH. Fluctuations of the belt α-helix are observed for closed conformations of the toxins and at acidic pH, while patches of more solvent-accessible residues appear under the same conditions in the core translocation domain HCNT. These findings suggest that, during translocation, the higher mobility of the belt could be transmitted to HCNT, leading to the favorable interaction of HCNT residues with the non-polar membrane environment.

"Water Association" Band in Saccharide Amorphous Matrices: Role of Residual Water on Bioprotection.

Saccharides protect biostructures against adverse environmental conditions mainly by preventing large scale motions leading to unfolding. The efficiency of this molecular mechanism, which is higher in trehalose with respect to other sugars, strongly depends on hydration and sugar/protein ratio. Here we report an Infrared Spectroscopy study on dry amorphous matrices of the disaccharides trehalose, maltose, sucrose and lactose, and the trisaccharide raffinose. Samples with and without embedded protein (Myoglobin) are investigated at different sugar/protein ratios, and compared. To inspect matrix properties we analyse the Water Association Band (WAB), and carefully decompose it into sub-bands, since their relative population has been shown to effectively probe water structure and dynamics in different matrices. In this work the analysis is extended to investigate the structure of protein-sugar-water samples, for the first time. Results show that several classes of water molecules can be identified in the protein and sugar environment and that their relative population is dependent on the type of sugar and, most important, on the sugar/protein ratio. This gives relevant information on how the molecular interplay between residual waters, sugar and protein molecules affect the biopreserving properties of saccharides matrices.

Thermodynamics and Kinetics of Ion Permeation in Wild-Type and Mutated Open Active Conformation of the Human α7 Nicotinic Receptor.

Molecular studies of human pentameric ligand-gated ion channels (LGICs) expressed in neurons and at neuromuscular junctions are of utmost importance in the development of therapeutic strategies for neurological disorders. We focus here on the nicotinic acetylcholine receptor nAChR-α7, a homopentameric channel widely expressed in the human brain, with a proven role in a wide spectrum of disorders including schizophrenia and Alzheimer's disease. By exploiting an all-atom structural model of the full (transmembrane and extracellular) protein in the open, agonist-bound conformation we recently developed, we evaluate the free energy and the mean first passage time of single-ion permeation using molecular dynamics simulations and the milestoning method with Voronoi tessellation. The results for the wild-type channel provide the first available mapping of the potential of mean force in the full-length α7 nAChR, reveal its expected cationic nature, and are in good agreement with simulation data for other channels of the LGIC family and with experimental data on nAChRs. We then investigate the role of a specific mutation directly related to ion selectivity in LGICs, the E-1' → A-1' substitution at the cytoplasmatic selectivity filter. We find that the mutation strongly affects sodium and chloride permeation in opposite directions, leading to a complete inversion of selectivity, at variance with the limited experimental results available that classify this mutant as cationic. We thus provide structural determinants for the observed cationic-to-anionic inversion, revealing a key role of the protonation state of residue rings far from the mutation, in the proximity of the hydrophobic channel gate.

Hepatitis B protein HBx binds the DLEU2 lncRNA to sustain cccDNA and host cancer-related gene transcription.

OBJECTIVE: The HBV HBx regulatory protein is required for transcription from the covalently closed circular DNA (cccDNA) minichromosome and affects the epigenetic control of both viral and host cellular chromatin. DESIGN: We explored, in relevant cellular models of HBV replication, the functional consequences of HBx interaction with DLEU2, a long non-coding RNA (lncRNA) expressed in the liver and increased in human hepatocellular carcinoma (HCC), in the regulation of host target genes and the HBV cccDNA. RESULTS: We show that HBx binds the promoter region, enhances the transcription and induces the accumulation of DLEU2 in infected hepatocytes. We found that nuclear DLEU2 directly binds HBx and the histone methyltransferase enhancer of zeste homolog 2 (EZH2), the catalytic active subunit of the polycomb repressor complex 2 (PRC2) complex. Computational modelling and biochemical evidence suggest that HBx and EZH2 share two preferential binding sites in DLEU2 intron 1. HBx and DLEU2 co-recruitment on the cccDNA displaces EZH2 from the viral chromatin to boost transcription and viral replication. DLEU2-HBx association with target host promoters relieves EZH2 repression and leads to the transcriptional activation of a subset of EZH2/PRC2 target genes in HBV-infected cells and HBV-related HCCs. CONCLUSIONS: Our results highlight the ability of HBx to bind RNA to impact on the epigenetic control of both viral cccDNA and host genes and provide a new key to understand the role of DLEU2 and EZH2 overexpression in HBV-related HCCs and HBx contribution to hepatocytes transformation.

Understanding and Controlling Food Protein Structure and Function in Foods: Perspectives from Experiments and Computer Simulations.

The structure and interactions of proteins play a critical role in determining the quality attributes of many foods, beverages, and pharmaceutical products. Incorporating a multiscale understanding of the structure-function relationships of proteins can provide greater insight into, and control of, the relevant processes at play. Combining data from experimental measurements, human sensory panels, and computer simulations through machine learning allows the construction of statistical models relating nanoscale properties of proteins to the physicochemical properties, physiological outcomes, and tastes of foods. This review highlights several examples of advanced computer simulations at molecular, mesoscale, and multiscale levels that shed light on the mechanisms at play in foods, thereby facilitating their control. It includes a practical simulation toolbox for those new to in silico modeling.

Closed-Locked and Apo-Resting State Structures of the Human α7 Nicotinic Receptor: A Computational Study.

Nicotinic acetylcholine receptors, belonging to the Cys-loop superfamily of ligand-gated ion channels (LGICs), are membrane proteins present in neurons and at neuromuscular junctions. They are responsible for signal transmission, and their function is regulated by neurotransmitters, agonists, and antagonists drugs. A detailed knowledge of their conformational transition in response to ligand binding is critical to understanding the basis of ligand-receptor interaction, in view of new pharmacological approaches to control receptor activity. However, the scarcity of experimentally derived structures of human channels makes this perspective extremely challenging. To contribute overcoming this issue, we have recently reported structural models for the open and the desensitized states of the human α7 nicotinic receptor. Here, we provide all-atom structural models of the same receptor in two different nonconductive states. The first structure, built via homology modeling and relaxed with extensive Molecular Dynamics simulations, represents the receptor bound to the natural antagonist α-conotoxin ImI. After comparison with available experimental data and computational models of other eukaryotic LGICs, we deem it consistent with the "closed-locked" state. The second model, obtained with simulations from the spontaneous relaxation of the open, agonist-bound α7 structure after ligand removal, recapitulates the characteristics of the apo-resting state of the receptor. These results add to our previous work on the active and desensitized state conformations, contributing to the structural characterization of the conformational landscape of the human α7 receptor and suggesting benchmarks to discriminate among conformations found in experiments or in simulations of LGICs. In particular key interactions at the interface between the extracellular domain and the transmembrane domain are identified, that could be critical to the α7 receptor function.

Bioprotection Can Be Tuned with a Proper Protein/Saccharide Ratio: The Case of Solid Amorphous Matrices.

Saccharides, and in particular trehalose, are well known for their high efficiency in protecting biostructures against adverse environmental conditions. The protein dynamics is known to be highly inhibited in a low-water trehalose host medium, the inhibition being markedly dependent on the amount of residual water. Besides hydration, the protein/sugar ratio is expected to affect the properties of saccharide amorphous matrices. In this work, we report an infrared spectroscopy study in dry amorphous matrices of various sugars (the disaccharides trehalose, maltose, sucrose, and lactose, and the trisaccharide raffinose) containing myoglobin, at different protein/sugar ratios. We analyze the stretching band of the bound CO molecule and the water association band. Such bands have already been successfully exploited for the simultaneous study of thermal evolution of a matrix and embedded protein. The results show a high dependence of protein and matrix signals on the protein/sugar ratio, the system behavior evolving from situations where (i) the protein slaves the matrix to (ii) protein ↔ matrix coupling/uncoupling, then to (iii) the matrix slaving the protein, with increasing sugar concentration. This supports a mutual protein ↔ matrix structural and dynamic influence in low hydrated systems, indicating that the protein/solvent master and slave paradigm does not strictly hold, but the mutual relationship depends on the relative concentrations. Furthermore, for each sugar, an optimal protein/sugar concentration ratio can be identified, which maximizes the protein preservation; under such a condition, the water content is minimal.

Biopreservation of Myoglobin in Crowded Environment: A Comparison between Gelatin and Trehalose Matrixes.

Biopreservation by sugar and/or polymeric matrixes is a thoroughly studied research topic with wide technological relevance. Ternary amorphous systems containing both saccharides and proteins are extensively exploited to model the in vivo biopreservation process. With the aim of disentangling the effect of saccharides and polypeptidic crowders (such as gelatin) on the preservation of a model protein, we present here a combined differential scanning calorimetry and UV-vis spectrophotometry study on samples of myoglobin embedded in amorphous gelatin and trehalose + gelatin matrixes at different hydrations, and compare them with amorphous myoglobin-only and myoglobin-trehalose samples. The results point out the different effects of gelatin, which acts mainly as a crowding agent, and trehalose, which acts mainly by direct interaction. Gelatin is able to improve effectively the protein thermal stability at very low hydration; however, it has small effects at medium to high hydration. Consistently, gelatin appears to be more effective than trehalose against massive denaturation in the long time range, while the mixed trehalose + collagen matrix is most effective in preserving protein functionality, outdoing both gelatin-only and trehalose-only matrixes.

The water association band as a marker of hydrogen bonds in trehalose amorphous matrices.

The relevant role played by residual water in modulating the dynamics and structure of a protein, a matrix and their coupling has been thoroughly studied in bioprotective amorphous saccharide matrices via experiments and simulations. In order to better characterize this residual water and the hydrogen bond structures in which it is involved, in this work infrared spectroscopy experiments are conducted on trehalose-water systems. The properties of water are inferred from the study of a peculiar infrared band, the water association band, which we exploited as a marker of the hydrogen bonds in which water is involved. Our aim was the identification of populations of water molecules, which give rise to the different components to which the water association band can be easily decomposed. The attribution of these components to families of water molecules is accomplished by studying the band behaviour with a suitable use of Hofmeister salts, known to have a structure-making or structure-breaking activity, and therefore able to modify the hydrogen bond network by enhancing or depressing the local order. The results allow ascribing, in almost all samples, five band components to either a chaotropic or kosmotropic environment, and further define two of them as bulk-like or ice-like water. The characterization of different components enables the use of this band as a tool to deepen the knowledge of other low-water hydrated matrices with a new approach. A differential analysis of peak frequencies and populations of the components in a bulky system, containing or not embedded components or interfaces (e.g. proteins, polymers, surfaces or even massive cosolutes), makes it possible to draw information on the properties of hydrogen bonds which are formed in the investigated systems.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation.

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking. Homology modeling and Molecular Dynamics simulations are valuable tools to predict structures of unknown proteins, however, for the case of human nAchRs, they have been unsuccessful in providing a stable open structure so far. This is due to different problems with the template structures: on one side the homology with prokaryotic species is too low, while on the other the open eukaryotic GluCl proved itself unstable in several MD studies and collapsed to a dehydrated, non-conductive conformation, even when bound to an agonist. Aim of this work is to obtain, by a mixing of state-of-the-art homology and simulation techniques, a plausible prediction of the structure (still unknown) of the open state of human α7 nAChR complexed with epibatidine, from which it is possible to start structural and functional test studies. To prevent channel closure we employ a restraint that keeps the transmembrane pore open, and obtain in this way a stable, hydrated conformation. To further validate this conformation, we run four long, unbiased simulations starting from configurations chosen at random along the restrained trajectory. The channel remains stable and hydrated over the whole runs. This allows to assess the stability of the putative open conformation over a cumulative time of 1 µs, 800 ns of which are of unbiased simulation. Mostly based on the analysis of pore hydration and size, we suggest that the obtained structure has reasonable chances to be (at least one of the possible) structures of the channel in the open conformation.

Ultrafast myoglobin structural dynamics observed with an X-ray free-electron laser.

Light absorption can trigger biologically relevant protein conformational changes. The light-induced structural rearrangement at the level of a photoexcited chromophore is known to occur in the femtosecond timescale and is expected to propagate through the protein as a quake-like intramolecular motion. Here we report direct experimental evidence of such 'proteinquake' observed in myoglobin through femtosecond X-ray solution scattering measurements performed at the Linac Coherent Light Source X-ray free-electron laser. An ultrafast increase of myoglobin radius of gyration occurs within 1 picosecond and is followed by a delayed protein expansion. As the system approaches equilibrium it undergoes damped oscillations with a ~3.6-picosecond time period. Our results unambiguously show how initially localized chemical changes can propagate at the level of the global protein conformation in the picosecond timescale.

Conformational changes in acetylcholine binding protein investigated by temperature accelerated molecular dynamics.

Despite the large number of studies available on nicotinic acetylcholine receptors, a complete account of the mechanistic aspects of their gating transition in response to ligand binding still remains elusive. As a first step toward dissecting the transition mechanism by accelerated sampling techniques, we study the ligand-induced conformational changes of the acetylcholine binding protein (AChBP), a widely accepted model for the full receptor extracellular domain. Using unbiased Molecular Dynamics (MD) and Temperature Accelerated Molecular Dynamics (TAMD) simulations we investigate the AChBP transition between the apo and the agonist-bound state. In long standard MD simulations, both conformations of the native protein are stable, while the agonist-bound structure evolves toward the apo one if the orientation of few key sidechains in the orthosteric cavity is modified. Conversely, TAMD simulations initiated from the native conformations are able to produce the spontaneous transition. With respect to the modified conformations, TAMD accelerates the transition by at least a factor 10. The analysis of some specific residue-residue interactions points out that the transition mechanism is based on the disruption/formation of few key hydrogen bonds. Finally, while early events of ligand dissociation are observed already in standard MD, TAMD accelerates the ligand detachment and, at the highest TAMD effective temperature, it is able to produce a complete dissociation path in one AChBP subunit.

Multiphoton absorption of myoglobin-nitric oxide complex: relaxation by D-NEMD of a stationary state.

The photodissociation and geminate recombination of nitric oxide in myoglobin, under continuous illumination, is modeled computationally. The relaxation of the photon energy into the protein matrix is also considered in a single simulation scheme that mimics a complete experimental setup. The dynamic approach to non-equilibrium molecular dynamics is used, starting from a steady state, to compute its relaxation to equilibrium. Simulations are conducted for the native form of sperm whale myoglobin and for two other mutants, V68W and L29F, illustrating a fair diversity of spatial and temporal geminate recombination processes. Energy flow to the heme and immediate protein environment provide hints to allostery. In particular, a pathway of energy flow between the heme and the FG loop is illustrated. Although the simulations were conducted for myoglobin only, the thermal fluctuations of the FG corner are in agreement with the large structural shifts of FG during the allosteric transition of tetrameric hemoglobin.

Protein thermal denaturation and matrix glass transition in different protein-trehalose-water systems.

Biopreservation by saccharides is a widely studied issue due to its scientific and technological importance; in particular, ternary amorphous protein-saccharide-water systems are extensively exploited to model the characteristics of the in vivo biopreservation process. We present here a differential scanning calorimetry (DSC) study on amorphous trehalose-water systems with embedded different proteins (myoglobin, lysozyme, BSA, hemoglobin), which differ for charge, surface, and volume properties. In our study, the protein/trehalose molar ratio is kept constant at 1/40, while the water/sugar molar ratio is varied between 2 and 300; results are compared with those obtained for binary trehalose-water systems. DSC upscans offer the possibility of investigating, in the same measurement, the thermodynamic properties of the matrix (glass transition, T(g)) and the functional properties of the encapsulated protein (thermal denaturation, T(den)). At high-to-intermediate hydration, the presence of the proteins increases the glass transition temperature of the encapsulating matrix. The effect mainly depends on size properties, and it can be ascribed to confinement exerted by the protein on the trehalose-water solvent. Conversely, at low hydration, lower T(g) values are measured in the presence of proteins: the lack of water promotes sugar-protein interactions, thus weakening the confinement effect and softening the matrix with respect to the binary system. A parallel T(den) increase is also observed; remarkably, this stabilization can reach ∼70 K at low hydration, a finding potentially of high biotechnological relevance. A linear relationship between T(g) and T(den) is also observed, in line with previous results; this finding suggests that collective water-trehalose interactions, responsible for the glass transition, also influence the protein denaturation.

Myoglobin embedded in saccharide amorphous matrices: water-dependent domains evidenced by small angle X-ray scattering.

We report Small Angle X-ray Scattering (SAXS) measurements performed on samples of carboxy-myoglobin (MbCO) embedded in low-water trehalose glasses. Results showed that, in such samples, "low-protein" trehalose-water domains are present, surrounded by a protein-trehalose-water background; such finding is supported by Infrared Spectroscopy (FTIR) measurements. These domains, which do not appear in the absence of the protein and in analogous sucrose systems, preferentially incorporate the incoming water at the onset of rehydration, and disappear following large hydration. This observation suggests that, in organisms under anhydrobiosis, analogous domains could play a buffering role against the daily variations of the atmospheric moisture. The reported results are rationalized by assuming sizably different protein-matrix coupling in trehalose with respect to sucrose, analogous to the one suggested for the photosynthetic reaction centre from Rhodobacter sphaeroides (F. Francia et al., J. Am. Chem. Soc., 2008, 130, 10240-10246).

Mapping the network of pathways of CO diffusion in myoglobin.

The pathways of diffusion of a CO molecule inside a myoglobin protein and toward the solvent are investigated. Specifically, the three-dimensional potential of mean force (PMF or free energy) of the CO molecule position inside the protein is calculated by using the single-sweep method in concert with fully resolved atomistic simulations in explicit solvent. The results are interpreted under the assumption that the diffusion of the ligand can be modeled as a navigation on the PMF in which the ligand hops between the PMF local minima following the minimum free energy paths (MFEPs) with rates set by the free energy barriers that need to be crossed. Here, all the local minima of the PMF, the MFEPs, and the barriers along them are calculated. The positions of the local minima are in good agreement with all the known binding cavities inside the protein, which indicates that these cavities may indeed serve as dynamical traps inside the protein and thereby influence the binding process. In addition, the MFEPs connecting the local PMF minima show a complicated network of possible pathways of exit of the dissociated CO starting from the primary docking site, in which the histidine gate is the closest exit from the binding site for the ligand but it is not the only possible one.

Thermal denaturation of myoglobin in water--disaccharide matrixes: relation with the glass transition of the system.

Proteins embedded in glassy saccharide systems are protected against adverse environmental conditions [Crowe et al. Annu. Rev. Physiol. 1998, 60, 73-103]. To further characterize this process, we studied the relationship between the glass transition temperature of the protein-containing saccharide system (T(g)) and the temperature of thermal denaturation of the embedded protein (T(den)). To this end, we studied by differential scanning calorimetry the thermal denaturation of ferric myoglobin in water/disaccharide mixtures containing nonreducing (trehalose, sucrose) or reducing (maltose, lactose) disaccharides. All the samples studied are, at room temperature, liquid systems whose viscosity varies from very low to very large values, depending on the water content. At a high water/saccharide mole ratio, homogeneous glass formation does not occur; regions of glass form, whose T(g) does not vary by varying the saccharide content, and the disaccharide barely affects the myoglobin denaturation temperature. At a suitably low water/saccharide mole ratio, by lowering the temperature, the systems undergo transition to the glassy state whose T(g) is determined by the water content; the Gordon-Taylor relationship between T(g) and the water/disaccharide mole ratio is obeyed; and T(den) increases by decreasing the hydration regardless of the disaccharide, such effect being entropy-driven. The presence of the protein was found to lower the T(g). Furthermore, for nonreducing disaccharides, plots of T(den) vs T(g) give linear correlations, whereas for reducing disaccharides, data exhibit an erratic behavior below a critical water/disaccharide ratio. We ascribe this behavior to the likelihood that in the latter samples, proteins have undergone Maillard reaction before thermal denaturation.

Density functional theory study of the trans-trans-cis (TTC)-->trans-trans-trans (TTT) isomerization of a photochromic spiropyran merocyanine.

Density Functional Theory (DFT) calculations have been performed on the TTC-->TTT isomerization reaction of the open forms of the 1',3'-dihydro-8-bromo-6-nitro-1',3',3'-trimethylspiro[2H-1-benzopyran-2,2'-(2H)indole (8-Br-6-nitro-BIPS) system. The calculations were carried out in vacuo and in methylene chloride solution at different temperatures. Results are compared with the available experimental values of free energy difference and activation energy in solution.

GFP-mut2 proteins in trehalose-water matrixes: spatially heterogeneous protein-water-sugar structures.

We report investigations on the properties of nanoenvironments around single-GFP-mut2 proteins in trehalose-water matrixes. Single-GFPmut2 molecules embedded in thin trehalose-water films were characterized in terms of their fluorescence brightness, bleaching dynamics, excited state lifetime, and fluorescence polarization. For each property, sets of approximately 100-150 single molecules have been investigated as a function of trehalose content and hydration. Three distinct and interconverting families of proteins have been found which differ widely in terms of bleaching dynamics, brightness, and fluorescence polarization, whose relative populations sizably depend on sample hydration. The reported results evidence the simultaneous presence of different protein-trehalose-water nanostructures whose rigidity increases by lowering the sample hydration. Such spatial inhomogeneity is in line with the well-known heterogeneous dynamics in supercooled fluids and in nonsolid carbohydrate glasses and gives a pictorial representation of the sharp, sudden reorganization of the above structures after uptake <==>release of water molecules.